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I. 



OBJECTIVES 



The objectives of the project were: (1) to continue the 
visual, or graphical, analysis of the BRI0HA3 data (the 
initial effort of which as reported in Ref. 1) in an 
attempt to verify anomalous signal behavior already 
documented and to discovec any new anomalies whica might be 
present; (2) toward that end, to make whatever changes were 
necessary in both the programs and data handling procedures 
to meet either new system constraints or to enhance over-all 
project efficiency; and (3) to analyze the data from a 
statistical standpoint (via two entirely new computer 
programs) in an effort to determine how often and by what 
amount (in decibels) several different classes of signals 
change. 

The format of the data, particularly that of the 2.h 
minute time duration limitation (although there were several 
time-wise contiguous data sets), made long term analysis, 
both visual and statistical, impossible. In addition, the 
daily intercept times for the data sets precluded any 
significant effort belag made to correlate signal 
characteristics from even approximately the same time of the 
day . 
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II. 



INTRODUCTION 



It is well documented that as radio signals traverse the 
ionosphere (including signals bounced off of one of its 
layers) they ace often cnanged in a Banner that goes beyond 
simple attenuation due to the distance travelled. Both the 
make-up of the ionosphere itself (how many ions are present 
at any one time, the elusive Sporadic-E ? n enomenon , 
night-day asymmetry and the geomagnetic non-reciprocity to 
name but a few) and the more, dramatic disruptions of 
ionospheric behavior (solar flare caused magnetic storm 
activity, polar cap absorption and the auroral displays are 
three) team with the Bore common, every day paenomena of 
"Faraday Rotation" and noise-producing lightning storms to 
alter any signal with which any or all of taem come in 
contact. The results can os severe attenuation, pronounced 
amplification, "ducting", "whistlers", and fading (sometimes 
called the Dellinger 2£ fact) . References 2-5 were 
excellent sources for ionospheric phenomena information. 

It was felt that the application of graphical techniques 
to signals which had traversed the ionosphere aignt add some 
insight into the problem of exactly what happens to these 
signals. 

For the most part, the visual presentation techniques, 
shown in various stages in Fig. 1, used in this study were 
the same as those presented in Ref. 1. However, some 
changes have been made and will be discussed in detail in 
the appropriate section. as in Ref. 1, interesting sets of 
scans have been reproduced as graph plots (using the 
VERSATEC plotter) for further analysis and presentation in 
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igure 1 - ANOMALY DATA PROCESSING STAGES 
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this paper. All thirty-four of the previously unanalyzed 
complete data sets were studied with varying degrees of 
success concerning anomalies (or, least, perceived to be so) 
found therin. 

The swept-tuned receiver and associated equipment used 
to acquire and record the data sets is an operational 
receiver equipment used for otnec than experimental data 
gathering. Its proper operation requires frequent tuning of 
numerous sensitivity adjustments. How often and by what 
amount these attenuators and thresholds are changed have a 
profound effect on the performance of the equipment. 
Therefore, programs HIN/'dAX 1 and MIN/MAX 2 were written in 
an effort to accumulate enough statistical data from the 
sets to warrant the r ecommendation of some appropriate 
algorithm for controlling the attenuator adjustments. 
Section IV of this paper will address itself to that topic. 



A. NATURE OF THE DATA 



The data was gathered using a swept-tuned receiver with 
either an omnidirectional or 12 degree beam antenna. The 
data sets were produced oy sampling the receiver output at a 
25 Hz sampling rate as it scanned downward through an 390 
kHz band. Each data set consists of, at most, 3597 scans 
(each scan being composed of 318 contiguous 2.3 kHz wide 
bins moving down in frequency as the bin nuuoer increases) 
placed on 7194 cards. Bird one gives amplitude information, 
in 2 dB high quantization levels, for bins 1 through 160 
while information on bins 161 through 318 is found on card 
two. The quantization levels are given in 5-oit binary 
words with even numbered bin information found in rows 11 
through 3 while odd numoered bin information is in rows 5 
through 9. In addition, each card has the card number and 



set number punched in cd^ 12 for accounting purposes. With 
each set there is a' parameter card which delineates the date 
and time of signal acquisition, base frequency, ajitenna 
configuration and attenuation applied. 

Of approximately 150 data sets cecorded, forty-four were 
made available for this study. Figure 2 is a table of 
pertinent information concerning the forty-two complete sets 
analyzed in Ref. 1 and this paper. 

B. CHANGES IN THE GRAPHICAL ANALYSIS PROCEDURE 



Equipment reconf igur ations since the graphical analysis 
procedure was formulated and reported in Ref. 1 and efforts 
to make the data .handling more efficient demanded several 
changes both in procedure and in the two main analysis 
programs. Anomaly A and Anomaly B. Specific information 
concerning the workings of the two programs may ns found in 
Ref. 1, particularly in section II B. 

First of all, the sheer number of binary cards to be 

transcribed onto magnetic tape (nearly a quarter of a 

million) obviated the card-to-tape technique used 

previously. A much faster and far more reliaole procedure 

% 

was initiated wherein the IBM 360/57 system, belonging to 
the W. R. Church Computer Center, was used after normal 
working hours. The cards were first read onto 9-track tapes 
(800 BPI, even parity) and then transcribed to the 7-track 
tapes (556 EPI, odd parity) required for the XDS-9300 
computer. Partial dumps were made of all sets from both the 
9- and 7-track tapes to ensure that the cards read onto the 
tapes properly and that tae tape-to-tape translation was 
successful. The programs used to read the cards onto tapes, 
translate the 9- to 7-track and dump the tapes for 
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Figure 2 - DATA SET INFORMATION 
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quality control can be found at the end of this paper. 

This new system of data loading resulting in the 
requirement to change the sections of the Anomaly A and 
Anomaly B programs refering to data read-in. This was 
accomplished with little difficulty. 

Next, a metasymbol subprogram was added to each of the 
analysis programs to enable entire files (data sets) to be 
skipped when a specific file required was located other than 
first on a tape. By calling the routine "PORSCN" a 
keyboard-entered number of files (NPILE = ) would be passed 
over. This proved to be extremely valuaole as a time and 
labor saving device. The routine "BAK5CN" was also available 
for backing over entire files but chronic tape drive 
malfunctions negated the use of this even more valuable 
tool. 

One final major program change was dictated by a 
computer center equipment reconfiguration wherein the 
CALCOMP 563 plotter was no longer available. The new system 
to be used was the PDP-11 computer in conjunction with a 
VERSATBC plotter. This cnange required modifications of the 
analysis program sections dealing with the processing of the 
output data (desired plots from the AGT-1D screen), found at 
the very end of the main program, and the generation of a 
new subroutine , called the Header Subroutine, which packs 
the header data into the proper format for the plotting 
unit. In addition, a PDP-11 plot package (used to interface 
the data 4:ape with the plotter) had to be extensively 
modified to meet the demands of this anomaly data. Listings 
of the revised Anomaly A and B programs, the metasymbol 
program for calling "FORSCN" and "5AKSCN" and the PDP-11 
plot routine (called "aa.c 1 * for Anomaly A plots - this 
routine alone is listed since the one utilized for the 
Anomaly B graph plots differs only in a few of the 
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approximately 433 lines comprising "aa.c") can oe found at 
the back of this paper. Subroutines 3INP, 3INPUT , FNS and 
VCD were not changed and, hence, ace not listed in this 
paper but may be found ia Ref. 1. 

C. ANOMALY A AND 3 PL0T5 

The plots which have resulted from the graphical 
analysis portion of tais study require some prior 
explanation. Due to the data handling restrictions imposed 
by the XDS-9300 computer, only a small portion of the data 
in each set could be viewed at any one time. This resulted 
in some difficulty in visualizing what could be long-term 
(time wise) or wide (frequency wise) phenomena and demanded 
that each set be analyzed several times before the entire 
set had been viewed. 

The Anomaly A plots consist of twenty scans (signal 
wavefronts sampled 0.3'4 seconds apart - or further if scans 
were deliberately skipped to view long-term effects) with 
the most recent scan in time falling at the bottom of the 
picture. The scans could be brought onto the screen in any 
number up to, and including, twenty at a time. ^s shown in 
Fig. 3, the vertical axis is marked off in 5 dB increments 
which apply only to tae bottom scan (called ’’base scan") . 
The horizontal axis is divided into ten 23 kHz wide sections 
each composed of ten frequency bins 2.8 kHz wide for a total 
of 100 frequency bins and 230 kHz.* This width can be varied 
to other than 100 bins if desired. The leftmost oin (termed 
the "base bin") is selected by keyborad input, thus allowing 
the analyst to slide up and down the data set 
(frequency- wise) at will. The scale along the horizontal 
axis, then, consists of increments of ten bins starting with 
the "base bin" number at the origin. The term "oase freq" 
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refers to the base frequency of the data set and decreases 
in 2.8 kHz increments as the bii number increases. Shading 
has been provided on some of the plots to highlight the 
suggested anomaly. 



The Anomaly B plots are a combination of the Anomaly A 
plot technique and a modification which enables the analyst 
to view selected bins (signals) in the time vs. amplitude 
domain. At the bottom of each plot ace five scans (Anomaly 
A style) used to aid in the selection of the bins to be 
plotted in the upper traces (see Fig. 4). In addition to 
the restricted number of scans to be seen at one time (only 
5) , computer limitations required that the number of 
frequency bins presented be reduced to eignty. Other than 
those changes, discussion of the Anomaly A plots is also 
valid for the lower set of Anomaly 3 plots. The upper set 
of Anomaly 3 plots is a collection of five keyboard-selected 
signals cnosen from among the scans below. The traces are 
numbered one to five, from top to bottom respectively. The 
horizontal axis is a time line composed of eighty data 
points at least 0.04 seconds aparc (depending upon the 
number of scnas skipped, if any) with the oldest point at 
the extreme right side of the trace) with the hasn marks on 
the vertical axis refecing to 10 dB amplitude levels. 
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AMPLITUDE TIME (SCANS 0.04 SECONDS APART) 

(5 dB/DIV) 




Figure 3 - SAMPLE ANOMALY A PLOT 
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AMPLITUDE (10 dB/DIV) 




3 




(10 dB/DIV) 



Figure 4 - SAMPLE ANOMALY 3 PLOT 
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III. 3R APHICAL ANALfSIS 



A. DISCUSSION 



Prior to initiation of the graphical analysis process 
via programs Anomaly A and 3, an enumeration of many of the 
possible "anomalies" expected was made. This was not to 
imply either that signals which had oeea markedly affected 
by the ionosphere would only exhibit these graphical 
manifestations of anomalous performance or that all forms of 
signal anomalies would be seen by the analyst. Tnis 
enumeration procedure was merely an attempt to 'visualize" 
what effects those sigaal-changing ionosphecic phenomena 
discussed in the introduction might have on the data sets 
yet to be analyzed. Although no attempt was made to make 
this study's "anomalies" conform witn those identified in 
Ref. 1, there has been an attempt to verify any previously 
seen. 

The graphical effects on the signal waveforms expected 
were: 



1. a sudden enhancenent/depression of one or more 
signals (long or short term) ; 

2. an abrupt discontinuity across an entire 
wavefront or poction thereof; 

3. a signal "bum?" or "valley" moving either up 
or down in frequency as the wavefront 
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progresses (Fig. 4, Ref. 1 refers) ; 



4. a shift in frequency of one or more signals 
(either short term, long term or aperiodic in 
nature) ; 

5. an unexplainable enhancement of one or more 
signals with a concomitant depression of one 
or more signals at the same time in the same 
data set; 

6. a complete loss of continuity of the signal 
wavefront for variable periods of time; 

7. unexpected variations in one or more signal 
contours as the wavefront progresses; 

8. any out-of-the-ordinary graphical implication 
of anomalous performance. 

Effects 1, 2, 4, 5, and 8 were discerned in the present 
study while several of these, plus effects 3, 5, and 7, were 
reported in Ref. 1 . 

8. ANOMALY A ANOMALIES 



As shown in Ref. 1 and verified in this paper, most of 
the data sets (whan viewed on the graphics terminal) 
exhibited a great measure of uniformity in signal makeup in 
the time domain and did not yield any obvious graphical 
anomalies. In fact, fulLy fifty per cent of tna BRIGHAM 
data sets failed to provide any repoctaDle irregularities in 
signal pattern with the other fifty per cent giving fairly 
uniform signal wavefronts for most of the time recorded. 
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The latter group did, however, exhibit anomalous (or, at 
least, assumed to be so) characteristics which were 
verifications of performaace reported in Ref. 1 or were 
newly discovered. The uniformity mentioned above aided 
immeasureably in the detection of anomalies since even a 
slight wavefront variation could be discerned ia the midst 
of otherwise uniform signal progression. 

Figures 5 and 5 are two examples of minor 
discontinuities (affecting 30 to 40 bins - approximately 100 
kHz) which appear to be the result of a multi-signal 
amplitude enhancement for a very short period of time. This 
effect was found scattered through almost half of the 
thirty-four data sets viewed via Anomaly A. Figure 7 is an 
example of aperiodic discontinuities which ran across almost 
all of the entire data set bandwidth. This change did not 
occur over the entire time duration of the sat and, as 
mentioned above, could not be traced over its entire 
frequency spectrum. This type of discontinuity was also 
found in data sets 33 and 80. Both of the discontinuity 
examples were anticipated although tae latter form holds the 
most promise, anomaly-wise. 

Fading was another effect expected but tae examples 
found were somewhat surprising in that it was anticipated 
that the entire frequency bandwidth of a signaL wavefront 
(890 kHz) would be affected. This proved not to be the 
case, with all of the examples of fading implicating at most 
20 frequency bins at a time. Figure 3 shows the signal 
centered on bin 148 fading in and out, a process which 
lasted several seconds. Figure 9 shows a much smaller 
time-duration fading centered on bin 24. 

A third type of possioLe anomalous activity is that of 
the shifting, up or down in frequency, of either a single 
signal or entire 390 kHz wide signal scans. The precise 
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TIME 




Set 3, Base Scan 141, Base Freq 3.5 MHz, Base Bin 1 



Figure 5 - NARROW DISCONTINUITY EXAMPLE #1 
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TIME 




Set 3, Base Scan 1023, Base Freq 3.2 MHz, Base Bin 100 



Figure 6 - NARROW DISCONTINUITY EXAMPLE #2 
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Figure 7 - WIDE DISCONTINUITY EXAMPLE 
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Set 53, Base Scan 544, Base Freq 6.7 MHz, Base Bin 100 



Figure 8 - FADING EXAMPLE #1 
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TIME 




Set 33, Base Scan 270, Base Freq 5 MHz, Base Bin 1 



Figure 9 - FADING EXAMPLE #2 
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explanation of this phenomenon is not available but it mast 
be kept in mini that malfunctions in the receiver system 
itself could be the cause rather taan some 
ionosphere-r elated activity. Figure 10 shows a possible 
left shift (approximately 3 bins) and distortion of two 
side-by-side signals. Figures 11 and 12 are examples of 
apparent frequency shifts of 160 bins both of wnich turned 
out to be data set problems. Since each scan is composed of 
two data cards of binary coded amplitudes (bins 1 - 160 on 
card one and bins 161 - 313 on card two) , the order in which 
the two cards appear and, as importantly, the existence of 
both cards for each scan are absolutely vital to the correct 
representation of the signal wavefront. For example, if 
card one of a scan is missing, then card one of the next 
scan is taken as card two of the former. This has a domino 
effect which will graphically manifest itself as an apparent 
shift of 160 bins in all subsequent scans. Figure 11 
illustrates precisely tais effect. A manual search of the 
missing from the set. Figure 12 is an interesting twist on 
this problem ia that the signal apparently shifted 160 bins 
to the right (or, down ia freguency) and then, taree scans 
later, shifted right back again. Another search of the 
cards showed that ,as expected, two cards were missing, one 
to cause the initial 160 bin right shift and a second to 
return the data cards to their correct sequence three scans 
later. One last exampLe of signal shifting is by far the 
most interesting of the lot. As illustrated by Fig. 13, 
there is a shift down ia frequency of approximately 12 bins 
(34 kHz) affecting the entire 318 bins (890 kHz bandwidth) 
of apeciodically spaced scans. This phenomenon was seen 
over the entire length of the data set at apparently random 
intervals. Data sets 138 and 143 also exhibited this 
aperiodic frequency shifting, although tae numbec of bins by 
which the scans were shifted varied slightly over the three 
sets . 
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Figure 10 - SHIFTING EXAMPLE #1 
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Figure 11 - SHIFTING EXAMPLE #2 
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SHIFTING EXAMPLE #3 
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SHIFTING EXAMPLE #4 



One final category for possible signal anomalies 
consists of five different phenomena each of whica occured, 
with one exception, only once in the thirty-four data sets. 
Figure 14 is an example wherein three different 11 anomalies'* 
were present. (It should be noted that base scan = 
"various** for Figs. 14 and 15 means that each of the three 
sections of the waveforms in the two figures had different 
base scans and implies that this in no way detracts from the 
usefulness of the figures as examples of the affects in 
question.) Bins 1 througi 40 and 83 through 98 exhibit a 
very high degree of signal attenuation not seen in normal 
data sets (i.e. data sets not found to have either a 
receiver malfunction or test signal present) . Bins 41 
through 87 exhibit normal signal characteristics. Bins 101 
through 200 evidence a combination of unusual attenuation 
and normal signal activity, while bins 201 througi 270 show 
either a general distortion of the signals present or a 
fairly high amplitude level of ambient noise and/or 
distorted signals plus noise. None of these phenomena could 
be explained by receiver malfunction, data card sequencing 
or special signal characteristics (i.e. test signal, etc.). 
Figure 15 shows a data set exhibiting no continuity of 
signal whatsoever. None of the signal "bumps" can be 
followed in a uniform manner as time progresses. The theory 
that a serious receiver system malfunction was tie cause of 
the distortion was borne out by a call to the cognizant 
authority. This malfunction was also verified foe data sets 
25 and 54. The sawtooth effect seen at the middle, right 
hand side of Fig. 16 is a pattern not seen elsewiere either 
in this set or any of the other thirty-three sets analyzed. 
It shows wild fluctuations in seven adjacent amplitude 
points (2.8 kHz apart) in one scan, whereas the rest of the 
scan appears to be normal. No explanation for this 
phenomenon can be offered. Another unexpected anomaly was 
the two flat-topped signals, centered at bins 124 and 123, 
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Figure 14 - MULTIPLE ANOMALY EXAMPLE 
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Figure 15 - EQUIPMENT MALFUNCTION EXAMPLE 
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Set 3, Base Scan 1815, Base Freq 3.2 MHz, Base Bin 100 



Figure 16 - SAWTOOTH PATTERN EXAMPLE 
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found in Fig. 17. It appears to be the result of some 
special attenuation applied to those frequency bins alone, 
because, as can be seen by the rather "large'’ signal 
centered on bin 178, other signals appear not to be 
affected. The additional attenuation applied to this set 
during reception , 10 d3, is lower than levels applied to 
other sets which were not similarly affected. The best 
phenomenon observed in this final category is found in Fig. 
18. For some reason (other than the possioility of a 
different form of receiver malfunction than seen before) an 
otherwise normal data sat lost total signal continuity 
across the entire frequency Dandwidth witn subsequent scans 
exhibiting a highly distorted nature. While this could not 
have been caused oy a card sequencing problem, its source 
cannot be identified. 



C. ANOMALY B ANOMALIES 



The Anomaly B progran provided a different perspective 
from which to observe tie data sets since it presented two 
dissimilar plots on the screen at the same time. The lower 
traces consisted of five Anomaly A style three-dimensional 
plots, this time only eigity points (or bins) wide, while 
tne upper five traces were the two-dimensional time vs. 
amplitude plots of bins selected from the traces oelow. It 
was felt mat this different perspective would enhance the 
anomaly identification effort, that it did not was 
disappointing. 

Of the eighteen data sets analyzed (3, 7, 16, 18, 20, 
21, 33, 36, 37, 38, 40, 41, 51, 53, 51, 80, 125, and 139), 
Anomaly B proved useful in only two instances: (1) to 
provide a graphic presentation of the three signal types 
used in the Statistical Aaalysis Section of this paper; and 
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Figure 17 - ATTENUATED SIGNALS EXAMPLE 
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Figure 18 - COMPLETE DISCONTINUITY EXAMPLE 



(2) to illustrate one very interesting phenomenon which 
could not be seen by the Anomaly A presentation. As shown 
in Fig. 19 this phenomenon consists of an undulation in the 
amplitude of the five upper traces. This "lazy sinusoidal" 
variation appears to affect some signals more often than 
others. The term "lazy sinusoidal" should be used with care 
since it is not meant to imply that there is any strict 
periodicity to be found in the undulation or that the traces 
look exactly like sinusoids. However, many of the signals 
did evidence the more smooth (and, hence, sinusoidall 
variation seen in trace one in the figure. 

Although only relatively smooth waveforms are depicted 
in the figure, this effect was identified over a large range 
of signal types including those which keyed off and on; 
however, it was much easier to see this phenomenon on the 
smooth signals. This variation was seen in every data set 
analyzed although it occnred more freguently and noticeaoly 
in some than in others. 



The cause of this undulation could be: (1) variations in 
transmitter output power levels; (2) receiver system 
processing; or (3) ionosphere-related phenomena including 
Faraday Rotation, Sporadic-F and fading. 
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Figure 19 - AMPLITUDE UNDULATION EXAMPLE 
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iv. 



STATISTICAL ANALYSIS 



A. DISCUSSION 



In conjunction with the graphical analysis approach to 
identify signal anomalies attrioutable to ionospheric 
properties, an attempt was made to glean additional 
information from the BRIGHAM data by means of statistical 
algorithms; however, it should be noted taat no attempt was 
made to identify anomalies by this approach. Rather, the 
programs MIN/MAX 1 and MIS/MAX 2 were developed in an effort 
to determine how often and by what amount (amplitude in 
decibels) both signal wavefronts (scans) and individual 
signals (bins or sets of bins) change. It was hoped that 
this new data mignt prove useful in the ievelopement of an 
algorithm to adjust the myriad sensitivity controls found on 
the swept-tuned receiving equipment used to acguire and 
record the BRIGHAM data. Although this effort nas nothing 
directly to do with the signal analysis effort reported 
previously, it did take considerable time and effort to 
write the programs and create the desired data base. In 
addition, this evolution proved extremely helpful in tne 
overall effort of understanding signal characteristics, the 
primary goal of this study. 

For the purposes of tais inquiry, three different types 
of signals (one that stayed "oo" through most of the signal 
duration, one that stayed "off" through most of the signal 
duration and one that appeared to key on and off throughout 
most of the signal duration) were chosen from data sets 
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which were collected it various tines of the day. The 
quotation narks are used (hare and throughout the rast of 
the paper) to imply a general trend noted ovar the life of 
the signal not soma unvarying signal characteristic . Figure 
20 is an Anomaly B graph plot presenting classic examples of 
the three signal types. The "ON" signal (exhibiting that 
"lazy sinusoidal" amplitude variation mentioned in the 
graphical analysis sectioa) is trace 1, the "0N/3FF" signal 
is trace 2 and the "OFF" signal is trace 3. Traces 4 and 5 
are merely copies of trace 3 to fill out the plot. A table 
giving all pertinent data on the twenty-four signals 
analyzed may be found in Fig. 21. In the figure, set 
frequency refers to the base fregueacy of the set when the 
recording was initiate! while TDI means the time of 
interception for the data set. An attempt was made to 
analyze two signals from each of four different time 
categories (inreresting from the ionosphere's standpoint) 

, those being sunrise, nooa, sunset and midnight. With one 
exception, set 53 - 21D9LOCAL, this effort was ceaiized. 

B. STATISTICAL ANALYSIS PROGRAMS 

The first generation algorithm developed, called MIN/MAX 
1, provided for a scan-by-scan (down entire frequency 
bandwiaths - cross sections of many signals - at fixed 
times) analysis of a variable number of scans witiin a given 
set. MIN/MAX 1 identified each of the peaks and valleys 
along the scans, chose the highest peax and lowest valley 
and computed the average peak and average valley heights. 
The quantization levels given on the data cards were changed 
to decibels and the amount of attenuation applied at the 
time of acquisition was added back in to give received 
signal levels accurately. MIN/MAX 1 proved to be of only 
small value since signal wavefronts vice signals themselves 
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Figure 20 - THREE MIN/MAX SIGNAL TYPES 
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Figure 21 - MIN /MAX SIGNAL PARAMETERS 
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were analyzed. A program listing for MIN/MAX 1 
at the end of this paper. 



may be found 



The next generation effort, tabbed MIN/MAX 2, focused on 
specific signals (the three varieties delineated in the 
discussion section above) and aided a number of more 
valuable statistical features. Like its predecessor, 
MIN/MAX 2 identified each of the peaks and valleys, chose 
the highest peak and lowest valley and computed the average 
peak and average valley heights. 



Its next feature was the determination of two different 
types of signal changes: (1) the comparison of sample point 
amplitudes with a sliding refernce point to determine 
long-term changes in signal levels; and (2) the comparison 
of adjacent sample point amplitudes to determine short-term 
changes. The former was initiated with the bin 1 point as 
the reference point and all suosequent points compared with 
it until an amplitude was found which differed from the 
reference by a previously agreed upon number of decibels 
(for this test, 4, 6, 12 and 18 d3 changes were used in 
succession) . That point then became the reference point and 
the process continued for tne entire 2.4 minute length of 
the signal. The number of such long-term changes and the 
elapsed time between each change were recorded. This data 
could then be used to attain the goal of this effort, the 
gauging of how often and oy what amount the various signal 
types changed. 



A further refinement of 
of amplitude histograms o 
changes to illustrate the ef 
of decibel fluctuation wh 
both the frequency of cha 
elapsed times. Since the r 
the smallest elapsed time po 



the procedure was the creation 
f the elapsed times between the 
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eceiver sampling rate was 25 Hz, 
ssible between X dB fluctuations 
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was 0.04 seconds. A Dcief analysis of the first set of 
elapsed times lad to the division of the time line 
(abscissa) into twenty-one time bins with the inclusive 
boundaries given in Fig. 22 and labeled Si through N21. As 
a further refinement, oins N1, N2 and the first segment of 
N3 were divided into individual 0.04 second bins labeled 
SI-1 through N1-7, N2-1 througn N2-3 and N3-1 respectively. 
The two histograms gave an excellent perspective from which 
to watch the changes in signal characteristic as the amount 
of decibel fluctuation tested was stepped from 4 to 6 to 12 
and, finally, to 1 8 dB. 



The final output of MIS/MAX 2 was, in addition to the 
peak and valley statistics mentioned above, a listing of tne 
number of long-term changes that took place in eaca of the 
thirty-seven time bins of interest, the percentage of the 
number of changes under each time bin as compared with the 
total number of changes (this allowed for a more valid 
comparison of the over-all change picture when ai increase 
in decibel deviation reguired resulted in fewer changes 
tabulated in the bins) , the total number of long-term 
changes (called NtJM in tae output) and the totaL number of 
changes in adjacent sample points (the short-term changes 
mentioned above) . In addition, the mean and standard 
deviation for the elapsed times between long-term changes 
were computed. A listing for MIN /M AX 2 may re found at the 
end of this paper. 



C. OBSERVATIONS AND CONCLUSIONS 



There were a number of surprises which resulted from the 
initial computer output from MIN/MAX 2. First, the high 
numner of both long- and short-term changes identified was 
not expected. Out of the possible 3597 changes (there were 
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HISTOGRAM TIME BIN PARAMETERS 



at most that many scans, k sample points, available for 
each signal of interest) the number of changes noted in some 
cases reached as high as 2300 for the lower decibel 
fluctuation levels. Second, the elapsed time histograms 
tabulated showed a marked tendency to be heavily damped 
exponentials in nature with the 3.04 second bin (N 1 — 1 ) 
having the highest number of occurences throughout. And 
third, there was a high degree of overlap between the 
statistical charac teristics generated for the three signal 
types chosen indicating that there was no absolute way to 
categorize a signal as "OS", "ON/OFF" or "OFF" merely by 
looking at its statistics. What was expected (and did 
occur) was that, as tae number of decibels used for 
comparison increased, the exponential distributions "flowed" 
to the right indicating greater elapsed times between larger 
decibel fluctuations. 

Figures 23 - 29 represent the bulk of the meaningful 
data garnered from the HIN/HAX 2 application to the 
twenty-four signals of interest. In speaking of the 
characteristics which differentiate the three arbitrarily 
chosen and defined signal types. Figs. 23 - 28 provide those 
statistics which can be of value in coming to a logical 
conclusion concerning their differences. As would be 
expected from "ON/OFF" signals, they showed a greater 
difference between both average peak and valley values and 
between "maxmax" and "minnin" values than did the other two, 
implying more signal level fluctuation. The "3FF" signals 
consistently showed the smallest fluctuations as would be 
expected from a series of either highly decayed signals or 
high but fluctuating ambient noise.' Based on these 
statistics alone it wonld be nearly impossible to 
distinguish with a high degree of certainty between an "ON" 
and an "ON/OFF" signal given only the parameters in those 
figures. Figure 29 provides both a summary of the data 
concerning tne number of changes tabulated for each decibel 
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Figure 23 - "ON" SIGNAL PEAK-VALLEY STATISTICS 
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"ON/OFF" SIGNAL PEAK-VALLEY STATISTICS 
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Figure 25 - "OFF" SIGNAL PEAK-VALLEY STATISTICS 
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Figure 26 - "ON" SIGNAL ELAPSED TIME STATISTICS 
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Figure 28 - "OFF" SIGNAL ELAPSED TIME STATISTICS 
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Figure 29 - OVER-ALL dB VS. CHANGE DATA 
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comparison level and a general means of distinguishing 
between the three types of signals. Althouga there are 
definite differences between the average figures for the 
various dB levels and some differences in the ranges of 
values, the overlap encountered in the ranges mates signal 
differentiation by this means risky. In snort, what 
distinguished the "ON" from the "ON/OFF" and "OFF" signals 
was thoroughly subjective and only loosely corroborated by 
the statistics generated. 

Additional information can be retreived from Figs. 26 - 
28 which provide the number of changes (both long- and 
short-term) found when the decibel differential was "walked" 
from 4 to 18 dB and the percentage of long-term changes 
which fell into the most meaningful of the tairty-seven 
histogram bins (for elapsed times between long-term 
changes) . As is clearly shown, the number of cnanges, NOM 
and ADJ, decreased as the dB level increased while the 
percentage figures in columns N1 and N 1 — 1 evidenced a 
movement away from the first elapsed time bin. Both of 
these character istics imply that tae average time between 
signal level changes at the various dB settings increases as 
the number of dB increase. To be specific, by taking the 
data from Fig. 30 it could be calculated that tie average 
time elapsed between 4 and 6 dB fluctuations foe the tinree 
signal types were: 
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Figure 30 - CURVE-FITTING DATA 
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This data confirms taat tha "ON/OFF" signals fluctuated 
more rapidly than the other two categories while the "ON" 
signals proved to be more stable than the "OFF". 

These results imply that, if the agency that operates 
the svept-tuned receiving system used to acquire and record 
the BRIGHAM data desires to acquire FSK-type signals 
(roughly, the "GN' 1 variety) , it could get away with changing 
the sensitivity adjustments much less frequently than if 
Morse-type signals (rougaly, tha "GN/GF?" variety) were 
desired. This, plus the rough data provided in Fig. 23, 
could provide a good data base from which a dial controlling 
algorithm could be adjusted to meet the needs of the signal 
types of interest at any time. 

As was mentioned above, MIN/MAX 2 provided for the 
calculation of the mean and standard deviation of the 
elapsed times tabulated for eaca of forty-eight signals (the 
twenty-four signals of interest at both 4 and 6 dB 
comparison levels). Although the mean values wece accurate 
the standard deviation values were not since the formula 
used in their calculation is valid for only normal (or 
Gaussian) distributions. Since the elapsed time histograms 
were apparently exponential in nature, another form of 
variance determination was required. This effort called for 
a curve-fitting technique to identify the histograms as 
being of a particular distribution type, use of the 
chi-squared test for gooliess-of-f it and manipulation of the 
chosen distribution's mean and variance formulas to derive 
the information desired. Attempts were made to fit three 
distribution functions to the histograms (Gamma, Poisson and 
Geometric) with the Gamma function proving to be the most 
successful. The Gamma distribution function, with 
parameters 0<cX <“>=>, 0<^ <°° is defined as follows: 
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_ - 1 <X~1 -X/A 
(H *)j?) X 3 , X>3 



f(X) = 



0 , 



x<0 



£ 

Where, Mean = <*/3 , Variance . 

An algorithm was written which real in tha observed 
histogram frequencies, calculated the expected Gamma 
frequencies (using an International Mathematical Subroutine 
Library subroutine caLled MD3AM) and calculated the 
chi-squared value (essentially, the sum of the least squared 
differences between the oDserved and expected frequencies) 
for use in the goodness- of- fit test. The chi-squared value 
was then compared with tha appropriate value found in the 
chi-squared distribution table (appropriate as to degrees of 
freedom and significance level) ; if it fell below that table 
entry, the fit was considered good, if not the fit was bad. 
As discussed in a ef s. 7 and 8, a good rule of thumb is that 
the smaller the chi-squared value for a given number of 
degrees of freedom, the better the fit. As a word of 
caution, even an excellent fit (extremely small chi-squared 
value) does not guarantee that the distribution is really 
the one being tested, it merely states that the data do not 
present sufficient evidence to contradict the hypothesis 
that the histograms possess Gamma distributions. 



For the curve- fitting at hand, knowledge of tie mean and 
an iterative variation in c* permitted an iterative 
calculation of /S thus providing all of the necessary 
information for tie algocitnm discussed above. Tie process 
was refined until the smallest chi-squared value could be 
identified. It snould be noted that only the Ml through M15 
time bins were used since very few elapsed times ever fall 
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beyond five seconds and that, as the individual histogram 
composition dictated, either seven or three bins were 
actually used in the algorithm (following a suggestion that 
bins containing five or fewer counts be combined to enhance 
the accuracy of the test) . The bin limits implied that the 
degrees of freedom used were five or one respectively. 

Of the forty-eight histograms to be fitted, six had 
counts in too few bins to be useful (would have implied zero 
degrees of freedom - not found in the taole) , eigateen were 
oad fits and the remaining twenty-four provided good fits to 
the Gamma distribution. Figure 30 gives the results of the 
curve-fitting effort. Knowledge of £< and the mean implies 
knowledge of /3 which, in turn, implies (by definition) 
knowledge of the variance. a combination of taa mean and 
variance values should be useful in the determination of how 
often the sensitivity dials on the receiving system should 
be adjusted for the desired signals to be acquired. 

Figure 31 provides a breakdown of the C< values vs. 
signal type and amplitude change variable for the 
twenty-four histograms which "fit" the Gamma distribution. 
It was felt that tais would permit the grouping of 
histograms with respect to some combination of the three 
categories above. The six ranges shown in the figure imply 
that any grouping by £< alone would be tenuous indeed as 
they overlap quite heavily. The only discernaole trend is 
the apparent lowering of the average c<. figure as the 
amplitude change variable increases from 4 to 5 dB. This 
would seem to contradict the defining equations for the 
Gamma distribution whica imply that a lower <X implies a 
more "damped-exponential" form for the distribution (when, 
in fact, the more damped case is found at the 4 13 change 
variaole) . However, the shape of the distribution is also a 
function of /3 which, as defined in Sean = , tends to 
affect the distribution saape in a more radical fashion 
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Figure 31 - SIGNAL TYPE VS. ALPHA TABLE 
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than does 0^. . In fact, as tae change variable increases, 
the mean increases (sea Fig. 30) and c< seems to decrease 
(see Fig. 31) implying a very pronounced increase in p 
which, in the defining aguations, is found to be most 
prominent in the exponential teem; but, this implies a 
distribution shape less heavily dampled in the lower (4 dB( 
change variable level. This trend of decreasing values, 
then, was to be expected and appears to be the only 
meaningful product of the analysis of the Fig. 31 lata. 

As an aside on the suoject of signal continuity during 
the brief recorded time span of the waveforms (around 2.4 
minutes) , it should be noted that data set number 7 bins 31 
("ON/OFF " ) and 91 ("OFF") exhibited very unusual 
characteristics when compared with the other signals in 
their classes. As shown on Fig. 27, bin 81 stiLl had over 
1000 changes at 13 dB difference while the other "ON/OFF" 
signals had fallen off to much smaller numbers. In fact, 
that bin had two changes at 54 dB when ever-increasing 
numbers were applied. The "OFF" signal, bin 91, evidenced a 
similar degree of consistency as it had many fewer 4 dB 
changes than the other "OFF 11 signals wnile falling off much 
more guickly as the number of decibels increased (see Fig. 
28). In fact, this signal had only a few more than 1700 2 
dB changes implying exceptional amplitude continuity when 
compared with any of the twenty-three other signals 
analyzed. To give a more grapnic representation of this 
phenomenon. Figs. 32 and 33 compare the long-term signal 
changes found in data sets 3 and 7 with the latter proving 
to be the more consistent of the two, although even set 3 
has a good measure of uniformity throughout its 3600 scans. 
The base scans for each of the five scan sets ( whicn show 
the beginning, middle and end of each data set - time-wise) 
are provided at the right of the pictures. 
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Figure 32 - LONG-TERM SIGNAL CONTINUITY - SET 3 
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Figure 33 - LONG-TERM SIGNAL CONTINUITY - SET 7 
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V. 



0MMENT5 AND RECOMMENDATIONS 



Of the two signal anomaly identification techniques. 
Anomaly A proved to be the most fruitful. Narrow and wide 
wavefront discontinuities, fading, frequency sniffing of 
entire wavefronts, frequency shifting of individual signals 
and several miscellaneous varieties were seen. Two of the 
above (discontinuities and shifting) were reported on in 
Ref. 1 with the rest being unique to this study. Several 
examples of pseu do- ano ma 1 i es were included to perform the 
dual function of warning against the false assumption of 
anomalous performance simply because the plot seemed to 
imply it and to provide pictures with the narrative in such 
cases. Analysis with Anomaly B proved to be relatively 
unproductive as only the interesting and frequently seen 
phenomenon of signal amplitude undulation (reported in Ref. 
1 and expanded upon here) - oeing noted. The diverse and 
amazingly complex nature of the ionosphere and its 
interaction with its environment precludes tie precise 
assignment of a cause to every effect (i.e., anomaly). 
Suffice it to say that, except for those anomalies which can 
be explained away by equipment malfunction, nothing has been 
found in this effort whicn will simplify the model of the 
ionosphere as aa ever-changing medium in wnici literally 
anything can happen to a radio signal which traverses it. 

The statistical analysis work performed, from creation 
of the MIN/MAX algorithms to the tedious curve-fitting 
effort, resulted in a wealth of predominantly raw data 
wnich, with the mean and variance values, should provide a 
foundation upon whicn a receiver system sensitivity 
dial-controlling a Igor it nm could be refined. It was noted 
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in the curve-fitting section that it is perhaps the weakest 
link in the statistical study since the chi-sguared test 
guarantees neither that the elapsed time nistograms are 
Gamma distributions nor that there night exist other 
distributions which would "fit" the histogram data as well 
as the Gamma function did. Howevec, a statistically 
reasonable argument can be made that the histograms {at 
least half of them) are pcobaoly Gamma in nature and, hence, 
the variance values calculated can be accepted as accurate. 

Should additional work be done with the BRIGHAM data, 
the following remark-s/reconnendations could prove useful: 

(1) all of the complete data sets in hand have been 
thoroughly analyzed witn program Anomaly A - additional 
effort along this line migat be a waste of time and energy; 

(2) althouga nearly twenty of the data sets in hand have 
not been analyzed with program Anomaly 3, the disappointing 
results thusfar might act as a warning as to expected future 
rewards from additional effort; 

(3) if Anomaly A and/or Anomaly B are to De used again, 
the 5AKSCN feature found in the METASYMBDL subruutine (it 
provides for tne backing over of a specified number of 
files, or data sets) should be made to work as a time/labor 
saving device - the problem is in the magnetic tape unit, 
not the program; 

(4) the recommendation (made initially in Ref. 1) for 
the creation of a computer algoritam to searcn out the 
anomalies in an automatic and, hence, more rapid manner than 
the visual approach is guardedly endorsed herein - although 
the speed at wuich the data sets could be analyzed would be 
improved by at least an order of magnitude, tne extreme 
complexity of trying to tell the computer what to look for 
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as wall as the loss of that very intangible but valuable 
human discretionary power could combine to made this effort 
a quagmire; 

(5) should additional work, be done utilizing tha MIN/MAX 
routines. Anomaly B should oe used first to enhance tha 
accuracy in salacting the different signal types to oe 
analyzed - that this was not done in the present study may 
have resulted in tha rather large ranges found in the data; 

(6) tha accuracy of any future curve-fitting work could 
be enhanced by the inclusion of the tfeibull and Lognormal 
distributions (or any of the level- crossing distributions) 
to the list of possible curves to which the elapsed time 
histograms are matched because sets of such data can be 
fitted to more than one distribution with equal accuracy and 
validity - this was not done here oecause tha suggestad 
distrioutions were complex eaougn and the time remaining 
short enough that a good effort could not have oeen put 
forth. 
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GLOSSARY 



1 . 



S^oca dic^E: thin horizontal layer (or 

electron density embedded in the 
resulting in signal deflection and/or 



patchi of high 
regular E layer 
absorption. 



2. Night^Day Asymmetry: Less attenuation at nignt than day 

due to difference in the ionizing effect of the sun 
resulting in discontinuities daring the sane ise/sunset 
periods . 



3* Geomagnetic iiEElFeciorocity: different attenuation 

effects on signals travelling west to east rather than 
east to west due to the earth's magnetic field lines. 



4. E2i§.E Cap i2§2.EEEi2! A|_ : absorption of radio signals 

passing the polar regions - due to a combination of 
solar flare activity and the shape of the earth's 
magnetic field lines. 

5. AEE2E§.i display: brigatly colored northern latitude 

phenomena disruptive to communications - also as the 
result of a combination of solar flare activity and the 
earth's magnetic field lines. 



6. Faraday Rotation: any linearly polarized wave travelling 

in the direction of a magnetic filed results in its two 
circularly polarizel components travelling at different 
velocities and thus the plane of polarization will 
rotate with distance. 
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7. Ducting: if a radio wm comas into contact with a ragion 

of inhomogeneous refractive indices the wave could be 
trapped between two layers and guided, as in a leaky 
wave guide, away froa its intended destination. 

8- Fading gDa 1L inger Effect^: sudden ionospheric 

disturbance (S. I. D.) which produces a compLate radio 
•’fade out" lasting from a few minutes to an hour or 
more - caused by soLar flare activity. 



69 



LIST DF REFERENCES 



1. Stapleton, D. V., High Freyuen cy Propagation An o mal ies 

M.S.EE Thesis, U. S. Naval Postgraduate School, 

Monterey, 1973 

2. Alpert, I. L . , la gig Nave Propagation and the 

l 2 [L 2 §£here i 7o 1^ 1, p. 192 to 397, Consultants 

Bureau, 1974 . 

3. Budden, K. 3., Radio Waves in the I gngspnec ? , d. 1 to 
542, University Press, 1961. 

4. Serson, N. C. , Editor, Radi 2 Absorption 12 the 

l222§2ll§£§/ P* 1 to 29, 1 06 to 317, Pergaaon Press, 
1962. 

5. Jordan, E. C. and Balmain, K . 3., romagnet ic 

Waves and Radiating Systems , p. 557 to 699, 

Prentice-Hall, Inc., 1968. 

6. Ratcliffe , J. A., An Introduction to the I gno sphere 

and Mag ne tosphere, o. 1 to 131, University Press, 

1972. 

7. Breiman, L . , Statistics:. With a 7i eg To wards 

Applications, p. 175 to 217, Houghton Mifflin Co., 

1973. 

8. Freund, J. E. , H2iE£2 Elementary Statistics, p. 272 

to 293, Prentice-Hall, Inc., 1965. 



70 



PROGRAM TO READ ONE CATA SET (APPROX 7200 CARDS) ONTG A 9-TRACK 
- IN THIS CASE, THE DATA SET WILL CCCUPY LABEL (CR REGION) CNE 
TAPE KFS61I WITH TAPE DENSITY 800 EPI AND EVEN PARITY. 
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PROGRAM TO TRANSFER DATA FROM A 9-TRACK TAPE TC A 7-TRACK TAPE - 
TFIS CASE, DATA SETS 1 TFRCUGH 4 FROM 9-TKACK TAPE NPS611 WILL 
CCCUPY LABELS 1 THROUGH 4 ON 7-TRACK TAPE ROWECB WITH TAPE DENSITY 
£56 EFI AND CDD PARITY. 
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SIBRCUTINE HEADER (KTAFE, IWCTH,NSCAN) 
INTEGER PBUF 

CINENSICN PBUF (2, 130) , IHBUF(260 ) 

ECU VALENCE (PBUF, IhBUF) 
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METASYMBCL PROGRAM USED TG PASS OVER (IN EITHER FORWARD OR REVERSE 
C IRE CT ION ) ENTIRE FILES OR A SPECIFIED NUMBER CF RECORDS WITHIN ANY 
FILE . 
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CDCVT CONVERTS A WORD TO BCC 

ALLS NONE 

ALLED BY ASGN » FCRSCNt BAKSCN 
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PLOT ROUTINE "aa.c" 



#de fine 


XINT 


320 


#de f i ne 


YINT 


8000 


#de f i ne 


DELY 


3200 


#de f i ne 


DEL X 


16 


#de f i ne 


AMPY 


8000 


^define 


NLINE 


1 


U de f i ne 


PI 


3. 


^define 


NP T 


150 


int zero 0 ? 




i n t one 


1 S 




struct { 


i n t 


np t ; 




i nt 


1 i n e # 




i nt 


1 d x # 




i nt 


Idy? 




i nt 


cot r ; 




i nt 


nsp t ; 




i nt 


spot l ; 




i nt 


spot 2 ; 




i nt 


spot 3 ; 




i nt 


h r # m i n # 


sec ? 


i nt 


n f t # 




i nt 


l aq # 




i nt 


s r # 




i nt 


id (3) ; 




i nt 


i p ; 




i nt 


sc i ; 




i nt 


mo # day # 


y r ; 


) head# 


*h ; 





int Sbuf (300) ; 
struct data ( 
int x t 
int yd; 

) ; 

struct data *d; 
int dbuf (40001; 
int np ; 
int nip 8 # 



int idev#odev#sc1x#scly#biasx,tdev? 

ma i n ( argc # a rgv ) 
int * *arqv » 

{ 

int i # j # n ? 
char * c s ? 

i f ( (tdev = oppn("/dev/sDo"» 1 )) < 0)< 
print f ( " cannot oo^n s p p " ) ; 
e x i t ( ) ; 

) 

if((pdev = oppn ( "/dpv/ rvp” # 1 ) ) < 0 ) { 
prints (" cannot op^n rvn O'*); 
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e x i t ( ) ; 

) 

if((idev - open ( " /de v / rm t 6 M * 0 ) ) < 0 ) { 
prjntf ("cannot open rmt6 #"); 
e * i t ( ) ; 

) 

i f ( a rqc > 1 ) { 

C 3 = arov (U ; 
n = 0 ; 

while ( *cs >= '0* *cs < = ’9*) 

n = n * 10 ♦ *cstt - *0’; 

n = * 1 1 ; 

for ( i =0 ; i<n; i++) 

inn(idev, dbuf,P00); 
printf ("number of records skipped, 

> 

whi le((inp(idev,hbuf,P00))>0) ( 
h = hbu f ; 

J = 0? 

for(i =0; i< h -> line? i+t)( 

n = i np ( i de v , 8-dbu f ( j 1 , P00); 

j = ♦ ?*(h->npt ) I 

) 

sclx = 030; sc 1 y = 03; 
b i as x = 50 ; 

np = h->npt * h - > 1 i n e ; 
scale( ) ; 
p ) o t ( ) ; 

c vers (pdev , 020 ) ; 
st t y (pdev , &one ) » 

} 



^define 


NB Y T 


2b4 


^define 


NBLK 


0 


tf de f i ne 


NSL 


1250 


a de f i ne 


DRAW 





st rue t i p t < 




i nt 


y * 


C h a r 


* x ob ; 


Char 


y i n c ; 


c ha r 


x d i r ; 


i n t 


y f ? 


i n t 


exp; 


i n t 


* f 1 ink; 


i n t 


*b 1 ink; 


) itabl7*NBYT),*io, 



char ob (NB Y T 1 ; 



Xd#",n) 
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i nt *dp ( 3000) 



plot ( ) 

< 



) 



S t rue t i p t * s ; 
i n t i r j f i e r 
sort ()? 

f o r ( i = 0; i < NBYT; i + ♦ ) 
pb ( i ] = 0 ; 

f or ( i = 0; i < NBLK? iff) 
w r i te(pdevrobfMBYT) ; 
ip r 0; i a r i t a b ; 

j = o; 

for(i - NSL? i > -l; i - - M 
whi1e(j < np R& *dp(J) 
i e = s i p ( dp ( j ) ) ? 
i f ( i e = = 1 ) 
re t urn ; 



//clear olot buffer 



//move to top of plot area 
//plot scan line 



= : i)i 

//set up plot point for interpolation 



J ♦ 

} 

n i b ( i ) J //set up plot buffer 

w r i t e ( pde v » pb » NR Y T ) ; / /p 1 o t line 



> 



i n t j s o r t ; 

so r t ( ) 

< 

reai ster i # k $ t ? 

i =0 ; 

f or (d = dbuf? d < &dbuf(np*2); d ♦ ♦ ) 
dp ( i ♦ ♦ 1 - & ( d -> yd); 

k -np ? 

while ( k =>> 1 ) < 

J s o r t + ♦ ; 

while ( j sor t ) { 
j sor t - 0 ; 

for ( i =0; i< (np-k ) ; i + O 

if ( * dp l i ) < *dp l i ♦ V ) ) { 

t = do f i 1 t 
dp ( i 1 =dp ( i ♦ k ) ; 
dp M +M =t ; 
j s o r t ♦ ♦ t 



} 
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s i p ( do i ) 

int * dp i ; 

< 

int x i , y i t x 1 , y 1 »xr,yr, incy»fyfdirx,cxf *Sf Icxr inc; 
i n t i f i e ; 
s-dpi ; 

1 c x - 1? 

yi = *dpi--; xi ” *doi; 



if(M S < Mdbuf (2*npl ) ) < 
x r = * 9 f + ; 

y r = ' * s ; 
i f (yr & DRAW) { 

i f ( ( i nc - y i - yr) >= 0 ) { 
if ( ( i nc y = i nc ) - = 0) 
di rx = nip; 
el sei 

for(i = 0; ( ( incy = ( inc/ (nip >> i))) 

di rx = (I << i); 

> 

f y - y r ; 

cx = 0200 ; 

lex = o; 

i e = stac*c(vi-incyfincyffyfxifdirx,cx); 
i f ( i == 1) 

returnt 1 ) ; 



= = 0 ) ; ilf 



) 



) 



} 

i f ( --dp i > = dbu f ) < 
i f (yi & DRAW) < 

y 1 - * d p i - - ; 

x 1 = *dp i ? 

if((inc = yi - yl) >= 0 ) { 
i f ( ( i nc y - inc) == 0 ) 

di rx = -nip; 

el sei 

for(i =0; ((incy - (inc/(nip >> i))) = = 0); MO 

dirx = -(1 << i); 

) 

f v = y 1 ; ' 
c x - lex; 

ie = stacHyi-incy» incy» fy» xi-1 ,dirx,cx); 
i f ( ie == 1 ) 

return( 1 ) ; 



> 



) 



) 

ret urn ( 0 ) ; 



stacMa,b»Cfdd,e, f ) 

int axbfCfddf^ff; 
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< 

i n t * 3 ; 
i n t i ? 

s t rue t i o t * 2 ? 

i a - > y - a ? 
i a - > vine = b ? 
i a - > y f = c ? 
i f (dd >= NBYT U dd < OH 

print f( "bad scale x=%dfl">dd)? 
e x i t ( ) ? 

) 

i a -> x pb - dd + pb ? 
ia -> xdfr - e ? 
ia - > exp = f ? 

i f ( ip == 0) ( 
i p= i t ab ? 

f p-> f 1 i nk = i o->b 1 i nk = 0 ? 
i a ♦ ♦ ; 

i a - > b 1 i n k = i p ; 
i a - > f 1 i n k = 0 ? 

} 

e 1 se i 

S = i a - > b 1 ink; 

s - > f 1 i n !r = i a i 
i f ( ia->f 1 ink == H){ 
s = ia? 

if ( ♦+ i a >= ^ i tab C 7 *NB Y T J ) < 

print f ("i tab overflow u " ) ? 
re t urn ( 1 ) ? 

) 

i a - > b 1 i n k = s ? 
i a- > f 1 ink=0? 

) 

else! 

s = ia -> flink? 
s - > blink = ia? 
ia - > flink = 0? 
i a = s ? 

} 

> 

) 



n i b ( s 1 ) 

i n t si? 

( 

i n t i f ) t n ? 
i n t * s ? s = ip? 

wh i 1 e ( s ) ( 

if(s - > yf < 0 ) 

i f ( (s-free(s) ) == 0) 
return; 

* s -> Xpb =S -> exp? 
i f ( s l = - s -> v ) < 

i - (Cn - s - > xdir) > 



//setup plottinq buffer 
//delete point 



//line break 
0 ? n : -n ) ? 
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//left 



i f ( n < 0 ) 

f o r ( j = 0 ; j < i ; I 4 4 ) { 
i f ( s -> c x o == 0 ) < 

* s -> xnb = ♦ 1 ; 
s - > exp = 1 ; 

} 

*s -> xnb =! (s -> exp 

> 

e I se 

f or ( } = o ; j < i ; i ♦ ♦ ) ( 

*s->xoh = ! ( s * > c x p =>> 

) 

s - > y = - s->yinc? 



) 

i f C S ) < = s->yf) //end 

s - > y f = - 1 r 
s = s->f 1 ink; 

> 

) 



f r ee ( s ) 

i nt * s ? 

< 

i n t * t # 
i n t W 

s t rue t i p t * z ; 

* s -> xpb - 0 ? 
i f (s*>b1 ink = - 0 ) f 
ip = s - > f 1 i n k ; 
i p - > b 1 i n k - 0# 

t = i p ; 

) 

e 1 se { 

t - s - > b 1 i n k ? 
t*>f link = 3 • > f 1 ink# 
t - s - > f 1 ink; 
i f (t == 0) ( 

s -> flink = i a ; 
i a = 3 ; 
return(t)? 

> 

t -> blink = s -> blink; 

} 

S-> flink = ia; 
s->b1ink = ia->b1ink; 
i a = s ; 
ret urn ( t ) ; 

} 



sc a 1 e ( ) 

< 

s t rue t data * s ; 

int dmf i » jfC^inx»miny; 

i nt dxl20)»dy(20j; 



= < < i ) ; 



l ) ; 



of point 
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s - dbuf; minx = miny ~ 077777; 
d= s ; 

f 0 r ( i = 0; i < np? i + + M 

minx = ((c = s - > x) < minx? c I 

v mi ny - ( ( c = S •> yd) < miny? c i 

s ♦ ♦ ? 

> 

if (minx > 0 ) 

minx • 0 ? 
if(miny > 0) 

miny - 0 ; 
s = dbu f ; 

for(i • 0; i < np; i ♦ f ) < 
dm = s *> yd M; 

S - > x - - minx; 
s - > yd = - miny; 
s - > x - / 9 C 1 X ; 

s - > x =♦ biasx; 
i f <s->x >= NRYT ) < 

printf ("overflow pb # " ) ; 
e x i t ( ) ? 

) 

s -> yd = / scly; 
i f <s->yd > NSL) 

s->yd = NSL - 1 ; 
s -> yd = *. 0177 776; 
s - > yd - ! dm ; 
s ♦ ♦ ; 

) 

) 

i np ( i d f , bu f , nby t e ) 

int i d f t *bu f $ nh v t p ; 

< 

int i / t f n , c ; 

S t rue t ( 

char cl*c2/c3,cu; 

) c f t 1 20 0 ) $ * s ; 
s = c f ; 

if((n = read(idf,cf#nbyte)) > OH 
f or ( i = 0; i < nbyte/y# i*+)< 
c - s - > c 2 << 2 ; 
t = c < < i o ; 
t =5 s -> c 3 << 6; 
t = ! s - > c<J; 
s 4 + ; 

bu ft i) = t ; 

} 

) 

) 

char *cbp / * 1 oc $ t M 1 32! , cb I 1 0) ; 



con v ( va 1 ) 

int v a 1 ; 



minx); 

miny) 
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< 



\ n t a ; 

i f ( a = val/10) 
conv ( a ) J 

*c bp ♦ + = val X 10 + *0» ; 



C one (cl r c 2 * n ) 
char *c2; 
i n t n,*cl; 

( • 

i n t i , m ; 

for(i = 0; i < n; i ++) < 



i f ( i == 


0) 


m - 


(*c 1 * 07700) >> 6 


else 




m - 


*cl *. 00 7 7 ; 


i f ( m = = 


0 12) 


*c2 


= »0'; 


i f ( m > = 


01 U m <= Oil) 


* c 2 


= M * + m - l ; 


i f ( m > = 


021 U m <= 031) 


*c2 


= • A ' ♦ m - 021 ; 


i f (m > = 


04 1 KK m <= 051) 


*c? = ' 


' J 1 ♦ m « 0 4 1 ; 


i f (m > = 


062 m <= 071) 


* c 2 


= 'S' + m - 062; 



c2 + ♦ ; 

} 

re t urn ( i ) ; 



s k i p ( cn t ) 

I nt ent ; 

( 

i n t i ; 

for(i = 0; i < ent? i+t) 
write(pdev,ph,?)J 

> 



clr() 

( 

i n t i ; 

for(i = 0; i < 13?? t ♦ + ) 

t Mil = • ’ ; 

1 1 M 3 1 1 = ' n • ; 

} 



mo v (tOf f r om , n ) 



107 



{ 



) 



char *tOf *f rom ; 
i n t n ? 

i n t i ; 

for(i = 0; i < n; H + ) 
♦ton s * f rom+ ♦ ? 
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SIMS THE VALLEY VALUES FOR AVERAGING PURPOSES 
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FRCGEAR TC PERPGRM A SERIES OF STATISTICAL CALCULATICNS WITH THE 
"ANCF'ALY" DATA ON A SELECTED BIN BASIS. 
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